Exchange engineering in 3d chains adsorbed on Cu2N/Cu(001) 
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Covalent substrates can give rise to a variety of magnetic interaction mechanisms among adsorbed 
transition metal atoms building atomic nanostructures. We show this by calculating the ground 
state magnetic configuration of monoatomic So? chains deposited on a monolayer of CU2N grown on 
Cu(OOl) as a function of d filling and of adsorption sites of the one dimensional nanostructures. 

PACS numbers: 



I. INTRODUCTION 

The study of the interactions which underlie the mag- 
netic ground state of deposited nanostructures consti- 
tutes an area of enormous interest nowadays, both, in 
applied and in basic research. The systems of interest 
can be atomically manipulated and deposited on differ- 
ent substrates by means of scanning tunneling microcopy 
(STM). For future spintronic applications the search is, 
among others, towards adequate films able to decouple 
the deposited nanostructures from the underlying metal- 
lic substrate, diminishing as much as possible the influ- 
ence of the latter on the electronic and magnetic proper- 
ties of adsorbed nano-objects. Actually, isolation is never 
complete and, the decoupling layers themselves play, in 
principle, a role to be determined on the magnetic inter- 
actions of the nanostructures. 

Copper nitride overlayers grown on different Cu sur- 
faces have been recently used as decoupling systems, be- 
ing it possible to develop stable thin films even at the 
monolayer thickness, while preserving an ordered geo- 
metric structure. Actually, the number of systems which 
fulfill these conditions seems to be limited and, there- 
after, copper nitride monolayers, called from now on 
CU2N, are currently being investigated both theoretically 
as experimentally.^ ^ 

Focusing on the magnetic interactions within the 
nanostructures, Mn chains of different lengths have been 
deposited by Hirjibehedin et al^ on Cu2N/Cu(001) at 
different adsorption sites by doing manipulation with 
an STM. Performing inelastic electron tunneling spec- 
troscopy (JETS) the authors concluded that the inter- 
atomic Mn-Mn interactions are antiferromagnetic with 
an interaction strength which depends on the deposi- 
tion configuration. There have been different theoret- 
ical contributions addressing these interactions .1^^^ In 
particular, in previous works we have studied the ef- 
fect of the surrounding geometry on the magnetic inter- 
atomic interactions for Mn and Cr chains deposited on 
Cu2N/Cu(001).l^ 

Recently the group of R. Wiesendanger has proposed 
to make use of the RKKY interaction among Fe adatoms 
deposited on a non magnetic metallic susbstrate, namely, 
the (111) surface of a copper monocrystal, to fabri- 



cate spin-based logic devices In fact, this group has 
achieved it to realize and develop a model system able 
to perform logic operations. The tip of an STM was 
used to deposit chains of antiferromagnetically coupled 
Fe adatoms and, using spin polarized scanning-tunneling 
spectroscopy (SP-STS), the magnetic interaction among 
the spins of the deposited adatoms was measured. The 
designed model device was able to transmit and process 
certain desired magnetic information. The authors sug- 
gest that combinations of antiferromagnetic and ferro- 
magnetically coupled chains could be used to realize new 
model gates by properly tuning the interatomic distances 
of the adsorbed atoms. 



In this contribution we show that using covalent de- 
coupling layers, such as the CU2N monolayers grown on 
Cu(OOl), it is possible to obtain different intrachain mag- 
netic interactions for an adsorbed atomic chain with 
a given d filling. This can be achieved not by tuning the 
interatomic distance of the adsorbed nanostructures but, 
by changing the adsorption sites of the transition metal 
atoms, which could be part of an eventual spintronic de- 
vice. The richness provided by the molecular substrate 
lies in the fact that the interatomic interactions within 
the deposited chains can change dramatically not only 
as a funtion of d-filling but also depending on adsorption 
configuration, while keeping the distance among chain 
atoms fixed. 



We present the results of a systematic study of the be- 
haviour of the magnetic interactions within Cr, Mn, Fe 
and Co chains deposited on the above mentioned molec- 
ular system. Four different chain adsorption geometries 
are considered, geometries which involve different mag- 
netic mediation processes, with the aim of highlighting 
the effect of the local and crystal environment on the 
interacting orbitals of the chains. We find that the 
adsorption sites of the chain atoms on the decoupling 
layer determine the kind of magnetic interaction involved, 
while the filling of the d-orbitals determines the final 
magnetic ground state. 



2 



II. CALCULATION DETAILS 

We perform ab initio calculations based on density 
functional theory (DFT) using the full potential lin- 
earized augmented plane waves method (FP-LAPW), as 
implemented in the WIEN2k codef^^The generalized gra- 
dient approximation (GGA) for the exchange and corre- 
lation potential in the parametrization of Perdew, Burke 
and Ernzerhof^ and the augmented plane wave- local or- 
bital (apw-lo) basis are used. The cutoff parameter which 
gives the number of plane waves in the interstitial region 
is taken as RmtKmax = 7, where K^ax is the value of the 
largest reciprocal lattice vector used in the plane waves' 
expansion and Rmt is the smallest muffin tin radio used. 
The number of k points in the Brillouin zone is enough 
in each case to obtain the desired energy and charge pre- 
cision, namely 10~^ Ry and 10~^ e respectively. 

The nitrogen atoms of the pristine substrate occupy 
fourfold coordinated hollow sites on the outermost metal- 
lic copper layer, forming a c(2x2) local structure.!^ This 
structure forms islands on the Cu(OOl) surface and, as 
mentioned in the introduction, the chains are built atom 
by atom on one of these CU2N islands. We simulate an is- 
land by considering a CU2N monolayer on Cu(OOl). This 
system is modeled with a supercell composed by a five 
layers slab separated by a vacuum region of 13.8 A. The 
thickness of the vacuum region is found to be sufficiently 
large to avoid interactions among subsequent slabs. To 
exploit inversion symmetry and to prevent unphysical 
multipoles, each slab has one CU2N monolayer on each 
side. 

The surface lattice constant follows pseudomorphically 
the underlying Cu(OOl) substrate and we use the exper- 
imental lattice constant of copper. We consider that the 
adsorbed 3d atomic chains are infinite long and they are 
deposited on both sides of the slabs. The nearest neigh- 
bour 3d-3d distance along the chains is 3.61 A in all the 
cases studied. In these supercell calculations, chains are 
periodically arranged 7.2 A apart from each other as- 
suming that the interchain interactions are much smaller 
than the intrachain ones. 

The positions of all atoms in the supercell of the pro- 
posed systems are allowed to move until forces are smaller 
than 0.1 eV/A. 

III. CALCULATIONS AND RESULTS 
A. The different surface structures 

For Cr, Mn, Fe and Co chains, we consider four dif- 
ferent stable adsorption geometries. Fig. [l] (a) displays a 
schematic top view of the Cu2N/Cu(001) substrate, while 
the different geometric chain configurations are depicted 
in Fig.[l](b-e). 

In the first adsorption geometry, the 3d chains are ar- 
ranged directly on top of Cu atoms of the CU2N mono- 
layer, in such a way that one N atom of the nitride lies 



inbetween each pair of 3d atoms along the chain. There 
are four lateral copper atoms of the CU2N monolayer sur- 
rounding each chain atom. This is one of the experimen- 
tally reported situationJ^ and we can it si. SeeFig.jlJb). 

In the second geometry, s2, each atom of the chain 
is deposited, as in the case of sl^ on top of Cu atoms 
but now every chain atom is laterally surrounded by two 
nitrogens, one at each side in the direction perpendicular 
to the chain, as shown in Fig. [l](c). There are also four 
lateral copper atoms of the CU2N substrate surrounding 
each chain atom. 

In the third case, called the s3 arrangement, the 3d 
atoms are located in hollow sites of the CU2N layer, as it 
is shown in Fig. [ijd). There are four surface Cu atoms 
and four nitrogens around each 3d one. 

The last configuration considered, the s4, is one in 
which the 3d atoms sit on top of nitrogens of the substrate 
as depicted in Fig.[l](e). Each 3d atom is surrounded by 
four coppers of the underlying substrate. This configu- 
ration has also been experimentally reported.^ 

Due to the reduced coordination number and the co- 
valent nature of the underlying CU2N network, the local 
atomic environment should play, a priori, a relevant role 
on the final electronic and magnetic structure of the de- 
posited chains, being it important to properly optimize 
them. We relax the different systems by starting always 
from the ferromagnetic spin configuration of the chain 
arrangement. 

For the different 3d fillings the final relaxed positions 
within each configuration are very similar, the differences 
among them being at most of 0.1 A. In the si geometry 
the N atoms, which are inbetween two 3d ones, change 
significantly their vertical distance with respect to the 
CU2N layer, relaxing outwards and forming a nearly lin- 
ear diatomic chain together with the 3d atoms. These 
last ones lie slightly higher (0.1 A) than the N atoms, 
the bond length between them is around 1.8 A for all d 
fillings considered. This is in agreement with previous 
results.^ The Cu atoms just below the diatomic chain 
are second nearest neighbors of the 3d ones. The dis- 
tance between them lies in the range 2.3-2 A A. There 
are four Cu atoms on the surface layer which are third 
nearest neighbours to a 3d one at 2.7-2.8 A distance. 

When relaxing the s2 structure, the two lateral N 
atoms of each 3d one move towards the chains forming 
N-3d bonds of lengths lying in the range 1.8-1.9 A. The 
distance between the 3d transition metal atoms and the 
Cu atom directly beneath them ends up being, after re- 
laxation, of 2.3-2.4 A. There are four lateral copper atoms 
at 2.9 A from each 3d one. 

The effect of relaxing the atomic positions in the s3 
configuration is that each 3d atom ends up having, in 
the CU2N monolayer, four copper atoms at 2.4-2.6 A and 
four nitrogens at 2.5 A . The copper atom beneath each 
3d one, which belongs to the metallic Cu(OOl) substrate 
is at 3.4-3.5 A of the corresponding chain atom. 

In the case of the s4 structure the distance between 
each 3d chain atom and the nitrogen sitting directly be- 
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Figure 1: (Color online) Schematic top view of the different adsorption configurations of the 3d chains on Cu2N/Cu(001) 
considered, (a) Pristine CU2N monolayer, (b)-(e) correspond to the si, s2, s3 and s4 arrangements. TM, N and Cu atoms are 
depicted by green, white and red spheres respectively 



low, ranges from 1.7 A to 1.8 A depending on d filling. In 
this case the N atoms move away from the CU2N mono- 
layer towards the 3d chain and the four nearest copper 
atoms are at 3.0-3.1 A distance from each chain atom. 

In all the adsorption geometries considered, the N 
atoms closest to the 3d chains have a tendency to relax 
towards the transition metal chain atoms. 



B. Bonding and electronic structure. 

The electronic density distribution of the systems un- 
der study depends on the adsorption sites of the chains. 
It determines the type of bonding and, as it will be seen 
in the next section, also the nature of the magnetic in- 
teractions. 

Analyzing the electronic density distribution, we find 
that there are three different ways in which the con- 
sidered adsorbed 3d chains bind to the substrate: i) 
forming a diatomic covalent one dimensional nanostruc- 
ture "metalically" bonded to the substrate, ii) building 
a monoatomic chain which binds covalently to the sub- 
strate and, iii) building a monoatomic chain which binds 
metalically to the substrate. 

In Fig. [2] we present electronic charge density plots for 
the si system along planes: (a) containing the deposited 
chain and (b) perpendicular to it and going through a 
chain atom. The si structure shows ID covalent bond- 
ing among each 3d atoms and the N ones which build 
the diatomic chain. The binding to the substrate of this 
covalent chain is achieved through copper atoms which 
originally belonged to the pristine CU2N monolayer and 
that now sit below the 3d atoms of the chain. This bind- 
ing can be considered as "metallic-like" when compared 
with the charge distribution inbetween Cu atoms of a Cu 
substrate. The bonding in this case belongs to type i). 

In the s2 configuration, the 3d chains are covalently 
bonded to the substrate through the two nitrogen atoms 
lying lateral to each chain one. These nitrogen atoms re- 
main being covalently bonded to the original monolayer. 
This can be concluded from Fig. [3} where electronic den- 




Figure 2: (Color online) Charge density plot for a Mn chain 
deposited on Cu2N/Cu(001) in the si configuration, (a) along 
the chain, (b) perpendicular to the chain. 



sity plots for different planes are being shown. In Fig. [3] 
(a) the selected plane contains a 3d chain and is perpen- 
dicular to the surface. In Fig. [3] (b) the plane is perpen- 
dicular to the chain and shows that the N atoms on each 
side, lateral to a chain atom, form covalent bonds with 
it. We are then in the presence of a type ii) bonding. 

In the s3 structure the adsorbed 3d chain sits in a 
metallic- like channel of the surrounding CU2N monolayer, 
the channel being built by the copper atoms nearest to 
the chain (case iii). The bond is metallic like as it can 
be observed in Fig. [4] Fig. [4] (a) displays the electronic 
density along a plane which contains the 3d chain and 
is perpendicular to the surface and Fig. |4] (b) shows it 
in a plane perpendicular to the chain which contains one 
3d atom. The local environment of the transition metal 
atoms is similar to the one they have when adsorbed di- 
rectly on Cu(OOl). Just for comparison, we present in 
Fig. |5] the electronic density distribution along and per- 
pendicular to a 3d chain absorbed on holes of the Cu(OOl) 
surface. We call this last system s5 and it is clear that 
the environment of the chain in s3 and s5 is very similar. 
We have checked that there is no covalent bond between 
the chain atoms and their neighboring nitrogens and this 




( a) (b) 



Figure 3: (Color online) Charge density plots for a Mn chain 
on Cu2N/Cu(001) in the s2 configuration (a) along the chain, 
(b) perpendicular to the chain and containing a chain atom. 



strengthens the fact that the chain presents only metallic 
bonds with its surroundings. 
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Figure 4: (Color online) Charge density plot for a Mn chain 
on Cu2N/Cu(001) in the s3 configuration, (a) along a plane 
paralell to the chain, (b) perpendicular to the chain and going 
through a chain atom. 

In the s4 configuration the 3d chain is covalently 
bonded to the CU2N substrate through the nitrogen atom 
sitting below each chain atom. This can be observed in 
the electronic density plots of Fig. [6j Fig. [g] (a) presents 
a plane containing the 3d chain, which is perpendicular 
to the surface and Fig. [g] (b) shows a plane perpendicular 
to the chain and to the surface, which goes through one 
atom of the chain. The binding in this arrangement is 
clearly of type ii). 

All the electronic density plots shown in this work are 
for Mn chains, but no relevant differences are observed 
for the other d fillings. 

Summarizing in si, s2 and s4 the 3d atoms of the 
chains are covalently bonded to their nearest neighbor 
nitrogen atoms, but, in si the chain binds metalically to 
the substrate, while in s2 and s4 it does it covalently. On 
the other hand, in s3 the 3d chain is metalically bonded 
to the CU2N substrate. 




( a) (b) 

Figure 6: (Color online) Charge density plot for a Mn chain 
on Cu2N/Cu(001) in the s4 configuration, (a) along a plane 
paralell to the chain (b) perpendicular to the chain and going 
through one chain atom. 



C. Magnetic configurations 

As stated in the introduction, among the possible mag- 
netic solutions, we consider only collinear ferromagnetic 
(FM) and antiferromagnetic (AFM) ones for each atomic 
adsorption configuration. We are aware of the fact that in 
those cases in which the difference in energy between the 
FM and AFM magnetic order is small, other interactions 
could turn out being important and that the ground state 
might be non collinear, but, in this work we are mainly 
concerned with general trends. 

In order to analyse the evolution of the magnetic 
ground state as a function of d band filling and adsorption 
sites of the chains, we calculate the difference in energy, 
AE^ between the AFM and FM configurations for all the 
arrangements. If this difference is negative the magnetic 
ground state is AFM, if not, it is FM. The results ob- 
tained for si, s2, s3 and s4 appear in Fig. [7| As it can 
be seen si has a completely different magnetic behaviour 
than the other three configurations. In si the magnetic 
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interaction goes from FM to AFM with increasing d fill- 
ing, while the reverse is observed for s2 and 5^, both, 
presenting similar trends. As discussed in the previous 
section, in the s2 and s4 arrangements the 3d atoms 
are covalently bonded to neighboring nitrogens, but in 
si there is a nitrogen atom inbetween two 3d ones. The 
nitrogen atom in the diatomic chain gives rise to a su- 
perexchange like magnetic coupling among the 3d atoms. 
This kind of adsorption situation, which enables this ni- 
trogen mediated interaction, has already been treated in 
the literature for Mn chains and for several dimers.l^HI^ 

In s2 and s4 the atoms of the 3d chains are covalently 
bonded to the substrate through their neighboring nitro- 
gen atoms and in both cases the magnetic ground states 
go from AFM to FM as a function of d filling. There is 
a direct exchange interaction along the chains and an in- 
direct interaction through the substrate. The question is 
which is more significant. If one compares the results for 
AE obtained in the case of the s2 and s4 configurations 
with those corresponding to the unsupported chains at 
the same interatomic distance, s6, the same trends for 
the magnetic interaction are found as a function of d 
filling. The energy differences between ferro- and anti- 
ferromagnetic configurations are also of the same order 
of magnitude. The results for AE for the unsupported 
chains together with those corresponding to s2 and s4 
are presented in Fig[7| (c). These results show that the 
direct interaction among the 3d atoms is comparable in 
the three cases. This is also apparent from the compari- 
son of the electronic density plots of Fig [s] (a) and Fig |6] 
(a) with that of an unsupported chain at the same in- 
teratomic distance. The indirect interaction through the 
substrate is certainly affecting the energy difference for 
each particular filling, giving rise to a switching of the 
solution from FM to AFM in the case of Mn. In the 
case of Mn the direct interaction is small and competes 
with the other interactions at the considered interatomic 
distance, but it does not change the general trends. 

Finally, s3 presents a different behavior as a function 
of d filling when compared with the other configurations 
and the results for AE point towards the preponderance 
of another type of magnetic interaction. The difference 
in energy AE, in this case, evolves as if the chains were 
adsorbed directly on hollow sites of Cu(OOl), that is, the 
so called s5 configuration. This is clear from Fig [t] (b) 
where we compare AE for s3 and s5 and from which 
it is seen that the order of magnitude and the sign of 
AE is the same for the four d filling studied. In the 
s3 arrangement the interaction among atoms of the 3d 



chains seems to be predominantly RKKY-like, as it is the 
case when 3d atoms are directly adsorbed on Cu(lll).'^ 
It is interesting to compare our results with those ob- 
tained for the exchange coupling among Ti atoms in 
Ref. 8 when depositing a monolayer of titanium on the 
same Cu2N/Cu(001)surface. The coupling obtained by 
the authors along an N axis and along a hollow axis are to 
be compared with the exchange energies obtained by us 
in the si and s2 configurations of the deposited chains. 
Following the trends already discussed, for Ti chains in 
the si and s2 configurations we would have obtained also 
ferromagnetic and antiferromagnetic exchange couplings, 
respectively. 

We have repeated the calculations within the DFT+U 
framework^^ using /7 = 5 eV for Cr and Mn and U = 2 
for Fe and Co. The results obtained for AE do no change 
the general trends obtained. 



IV. CONCLUSIONS 

Summarizing, the covalently bonded nitride substrate 
CU2N offers an interesting template to taylor magnetic 
interactions among transition metal atoms by taking as 
variable the adsorption sites. We have shown this by de- 
positing 3d chains on four different adsorption sites with 
the same interatomic distance among 3d atoms and for 
four different d fillings. Due to the different type of bond- 
ing present, we obtain that the 3d atoms interact magnet- 
ically through completely different mechanisms depend- 
ing on adsorption geometry. For a given d filling, the 
interaction can be strongly AFM or FM for the same in- 
teratomic distance along the chain. This versatility could 
be very promising regarding the future possibility of mak- 
ing spin engineering at the nanometer scale. We expect a 
similar richness in the magnetic interactions, when using 
other covalently bonded substrates. 
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